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SUMMARY 


The Optimal Control Model (OCM) of human response serves as a mechanism 
for generating sample time histories of human tracking error in different 
anti-tank systems. The systems under study Include TOW (Tube-Launched 
Optically Guided System) , DRAGON (Shoulder Mounted) and ITV (Improved TOW 
Vehicle). The model-generated trajectories are compared with field-test 
data across several dimensions Including time-domain (temporal) statistics, 
frequency content and subjective comparisons on individual runs. 

MODELING APPROACH 

The objective of this work is to develop a computerized model for gen- 
erating human tracking error time histories in several different anti-tank 
systems. The systems Include those in common use by the US Army such as 
TOW, DRAGON and ITV. A fourth system - GLLD (Ground Launched Laser Designa- 
tor) is similar to TOW and will not be discussed here. Of these systems, 

TOW and DRAGON are basically command llne-of-sight (LOS) , whereas ITV is a 
rate command system. The model that is developed must produce accurate fac- 
similes of tracking error over a wide range of target trajectories, from 
crossing (straight-line) motion to maneuvering motion. The model must be 
causal in the sense that future target motions are unknown at the present 
time. 


For the systems and target passes considered here, target motion is 
restricted to the azimuth axis, l.e. the gunner and target vehicle are both 
at the same ground level. Tracking error in elevation arises solely from 
the human's Inherent motor and observation randomness. Although the model 
that we have developed treats both axes, we discuss primarily the results 
for the azimuth axis here. A more complete discussion and presentation of 
the results may be found in Ref. [1]. 

The Optimal Control Model 

The Optimal Control Model of human response is used as the mechanism 
for building the anti-tank tracking model. The OCM technology has been 
successfully applied in numerous contexts including pilot control, anti-air- 
craft artillery, etc. The mechanics of using the OCM to generate sample 
path time-histories (as opposed to statistical measures) is described in 
Ref. [2]. Our application follows this approach, with minor modifications 
>".0 account for the dual-axis nature of the tracking task. The pertinent 
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equations for Monte-Carlo/Sample Path Simulation using the OCM may be found 
In Ref . [ 2 ] . 

Application of the OCM requires specification, for each given system- 
display-manipulator dynamics, of 1) the operator's task objectives In terms 
of a set of cost functional weights, 2) the parameters that define the oper- 
ator's Inherent limitations, and 3) an "Internal" model of the target dyna- 
mics. With these Items specified optimal control and estimation theory Is 
used to obtain the human's feedback strategy, and generate closed-loop per- 
formance results. 

1. Task Objectives: For a basic tracking task, wherein the human attempts 

to keep the error e(t) small, we use a cost functional 

J(u) * E{e^(t) + Q. e^(t) + Q. u^(t)} (1) 

6 U 

The weighting Q. reflects a human's subjective weighting on error rate. 
It Is Indicative of strategy, style or technique and could be associated 
with the type of training on a given system. The weighting on error 
rate Induces a first-order lag that Is associated with the neuro-motor 
system dynamics. For each system we select to yield Tjj»neuro-motor 
time-constant ■ .1 sec. The value of Q. Is to be determined for each 
system, based on data comparisons. 

2. Human Operator Limitations: The primary human operator limitations 

modeled In the OCM are those associated with perceiving displayed quan- 
tities and executing Intended control motions. The observational sub- 
model In the OCM assumes that the human observes tracking error 

yj^(t) • e(t) and tracking error rate y 2 (t) ■ e(t). However, the human 
perceives a delayed and noisy replica of these signals via 

ypi(t) “ ■ ■^>1 + Vyj^Ct - T) 1 - 1,2 (2) 

The function F(*) represents a vlsual/indlfference threshold of value 
aj^ ■■ 0.6 mr T display gain on error and &2 “ on error rate. The 

time-delay T ■ .15 sec. Each observation noise Is white with covariance 

Vyi(^) “ ^yi E{yj(t)}/ ATTN (3) 

Where « .Oln (-20dB) ; the attention allocations are assumed split 
0.8 for azimuth vs 0.2 for elevation.^ The nominal parameter values 
associated with the observational submodel are assumed fixed for all 
systems and target types considered. 

The neuro-motor submodel for generating human corrective Inputs is 
given by 


u(t) + u(t) - u (t) + V (t) 
c u 


A more precise model would be to employ dynamic attention allocation. 
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The quantity Uj,(t) is the "commanded" control input that is generated 
from the Kalman Filter/Predictor/Gains cascade; is the neuro-motor 
time constant. The white motor-noise Vy(t) consists of an additive plus 
a ratioed component for each axis. The covariance of the motor-noise 
for azimuth and elevation axes, respectively, is 
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The crossfeeds p^ and p£^ model uncertainty/randomness in one axis 
resulting from manipulator motion in the other axis. For the systems/ 
targets studied, elevation commands u^g are small relative to azimuth 
commands u^^ (recall target motion is in azimuth only). Thus, Pgg and 
p^ are not readily obtained from the available data. Therefore, we 
have assumed 


AA 


EE 


and 


AE 


(5) 


The remaining quantities V^, V°g. P^. Pg^ are system/manipulator depen- 
dent. Their values must be elicited from model — data comparisons. 
Finally, Eq(4) shows that the motor-noise scales with commanded control 
input Ug(t). In some Instances, e.g. command LOS systems, it is more 
natural for motor-induced randomness to scale with commanded rate u^(t). 


3. Target Submodel: In the present application of the OCM to track ground 

targets, target velocity §-f(t) and acceleration 0f(t) are generally 
small. Thus, we use a simple internal model for target motion 


x(t) - Wj(t) (6) 

where x(t) is the human's internal representation of target velocity. 
The "driving noise" W(j(t) has covariance 

cov (Wj(t)] ■ 3 0^(t) + a 0^(t) (7) 

Note that the "truth" model is x(t) ■ 6^(t) . The values selected for a 
and 3 are 

3 - 10^, cx - .015^ (a) 


These values are constant across all systems and targets studied. 

Data-Model Comparison Procedures 

The OCM can be used to generate, for a given system and target trajec- 
tory, an ensemble of tracking error time histories Ej, ■ {ej(t) ; j • 1,...M}. 
These model-generated runs may be compared against an ensemble of equivalent 
data trials, Ej. Clearly, it is the statistics of these two ensembles that 
one would wish to compare via model-data validation tests. Several modes of 
comparison are possible, ' discussed below. For consistency we have found 
it useful to remove ti .aporal mean Pej from each run prior to analysis. 
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( 9 ) 


Thus, 


1 ” 

Oj<ti) - ej(tj) - ! U.J - jj ^ 


where i * 1,...N are the sampled values of ej(t). This procedure removes 
random variations in signal mean, DC bias offsets in field recording equip- 
ment, biases in human aiming point, etc. 


1 . 


Ensemble Analysis: The ensemble mean and standard deviation can be com- 

puted in the usual manner, 

1 ^ 

y (t) ■ Ensemble mean ■ ^ 52 e.(t) (10a) 

e ” j ti 3 


a (t) > Ensemble SD 
e' 


1 

M-1 



(e^(t) - M^(t)] 


2 

( 


(10b) 


Since the sample runs have been rendered zero (temporal) mean we would 
expect iJg(t) = 0 if the ensemble was stationary. 


2 . 


Temporal Analysis: If the ensemble E is stationary, ~ 

Og(t) z constant so that the Information content in the temporal ensem- 
ble is reduced to a single number, i.e. RMS tracking error. A more 
direct way to obtain average RMS tracking error is to compute the tem- 
poral statistic 



( 11 ) 


for each run, and form the composite, M-run average, via 


o 

e 




( 12 ) 


However, unlike the ensemble analysis, Og is meaningful only in the 
stationary case. Its computation in the non-statlonary case is possible, 
but of dubious interpretation. 


3. Frequency Domain Analysis: The RMS temporal metrics give an indication 

of total error power, they do not Indicate how th*s power is distributed 
over the frequency range, whether there exists resonances, etc. To 
obtain these later indicators of system response we compute, from the 
temporal ensemble E <■ (ej(t), J ■ 1,...,M}, a frequency domain ensemble 
of normalized error PSD, E* ■ (e*(ijo), j ■ 1,...,M} where 

N 


e?((i>) - 


ej 


E 

i-l 


e^(t^) exp[-J(i)t^J 


(13) 


Normalization greatly reduces the sensitivity to motor-noise and Inter- 
subject variability since it considers only relative power distribution 
over (J. Note that the PSD computations are strictly valid for a sta- 


- 382 - 



tlonary ensemble; their computation is possible for any ensemble, of 
course, but the interpretation in the non-stationary case is dubious. 

The ensemble E* of normalized error PSD can be averaged (in the 
same manner as E) via equations similar' to (10a - 10b) to yield the 
ensemble mean yg*(u]) and ensemble SD, Comparison of model and 

data PSD statistics is thus possible and provides another, interesting, 
facet for model validation. 


RESULTS 


In this section model-data comparisons are given for the three anti-tank 
systems considered. In each case it is necessary to provide a description of 
the system-manipulator dynamics, and values for the motor-noise parameters 
and error rate weighting Q.. 


TOW System 

The TOW system is a command LOS system consisting of a launch tube plus 
sight mounted on a viscous (rate) damped turret. Thus the torque supplied by 
the operator to point the sight varies with sight (i.e. control) rate. The 
dynamical model used for the TOW system is 

T(s) - ;j:^ ; T - 0.1 sec. (14) 

These dynamics are chosen for convenience^, and are viewed as representative 
of the manipulator (arm-viscous mount) characteristics. 

In our data-model analysis of the TOW system we found Q. = 0 gave best 
match. The motor-noise in the OCM is assumed to scale with coimanded (i.e. 
LOS) rate, and the pertinent noise parameters (obtained from matching RMS 
scores) are 


v° - [.05, .01], - .005. - .015 

The temporal RMS statistics, computed via Eqs (11)-(12), for three dif- 
ferent ensembles are given in Table 1. The TF and TS ensembles correspond 
to crossing targets at a range of 3Km with 6- * 5.47 and ■ .55 mr/sec, 
respectively. These ensembles are stationary. The TM ensemble corresfonds 
to a set of 19 maneuvering trials. In these cases the target was approaching 
the gunner following a serpentine path with ^ mr/sec, - .5 mr/secr 

peak values. Each target pass was somewhat different; this ensemble is not 
statlcnary . 

The model-uata comparisons shown in Table 1 are excellent for both 
azimuth and elevation axis tracking (The numbers in parentheses are the 
computed standard deviations in the RMS tracking errors.) Only the eleva- 
tion SD is not well-matched for the 'maneuvering trials. This di.5ciepancy is 
expected to be corrected by a dynamic attentlonal submodel. A comparison of 

^At least first-order dynar.lcs are required by the OCM. 
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model-data normalized PSD ensemble statistics is shovn in Figs. 1-2 for the 
IF ensemble. The results are in excellent agreement. A t-test performed 
polntwise was used to confirm the equality of the PSD means at the 95% con- 
fidence level. 

TABLE 1 - COMPARISON CF TEMPORAL TRACKING ERRORS, TOW SYSTEM 

Azimuth SD Elevation SD 

Oata Model Data Model 

TF 23 0.106(0.017) 0.119(0.018) 0.060(0.009) 0.052(0.009) 

TS 19 0.057(0.005) 0.054(0.004) 0.030(0.006) 0.031(0.005) 

TM 19 0.35 (0.10) 0.35 (0.07) 0.14 (0.02) 0.C5 (0.01) 

Individual runs produced by model and data can also be compared subjec- 
tively. Fig 3 is a comparison of a model and a data run from the TF ensem- 
ble. Fig 4 likewise is a comparison of model-vs-data trials r one of the 
maneuvering target passes. The "eyeball" similarity is quite impressive. 


DRAGON System 


The DRAGON is a shoulder-mounted system that consists of a launch tube 
plus sight. The front of the tube is pivoted on a support; the rear part of 
the tube rests on the operator's shoulder. Thus, as the operator (usually 
in a seated position) tracks a crossing target he must continuously move his 
shoulder by leaning his torso more and more to one side. There are no dyna- 
mics asscclated with the system per-se. The only dynamics are those associ- 
ated with the operator's torso — i.e. the control "manipulator". These dyna- 
mics are approximated as 


T(s) 


s/S +1 


2 

^ 8 V 

+ 2 i 

^ — V 



' n ' 

' n 


(15) 


where i) ■ 11 1 rad/sec, 6 • 3 + 1 and s * ‘IS for tilt in the azimuth axis. 

n — ■ — 


The motor-noise for the DRAGON system is assumed to scale with commanded 
angle/body tilt as randomness increases greatly if one is required to track 
vdtlle leaning to one side. The motor-noise parameters are 


v” • (8, .25), p - .0001, p_. - .000015 
u u 

The weighting on error rate = 0. 


The tracking error data for DRAGON consisted of r 40 passes of a crossing 
target at 1km with |6_| " 10 mr/sec. In approximately 1/2 of the runs the 
target moved from right to left (DR); in the other runs motion was from left 
to right (DL) . A comparison of the temporal RMS statistics of isodel-vs-data 
is shown in Table 2. The results ure excellent, but this is not a stationary 
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ensemble as the motor-noise covariance increases during the course of a run! 
Indeed, Figs 5-6 show the true nature of the tracking error ensemble for 
model and data. This is a more meaningful comparison than is temporal RMS 
error . 

TABLE 2 - COMPARISON OF TEMPORAL TRACKING ERRORS, DRAGON SYSTEM 


M 

Azimuth SD 

Elevation SD 

Data 

Model 

Data 

Model 

DR 21 
DL 22 

7.35(3.3) 
7.2 (2.6) 

6.82(0.9) 

6.82(0.9) 

2.78(1.5) 

2.35(1.4) 

2.73(0.35) 

2.73(0.35) 


Comparisons of normalized error PSP for model and data provides another 
yardstick for judging the effectiveness of the OCM application. As noted 
earlier, interpretation of these results must he made cautiously as the en- 
semble is non-stationary. Nevertheless, we can consider this as the "aver- 
age frequency content" in the error waveforms. Figs 7-8 contain the model- 
data frequency comparisons for the azimuth axes. The results are excellent. 
A final model-data validation test is via the subjective comparison of indi- 
vidual tracking error time histories. Fig 9 shows a typical data run vs. a 
typical sample path from the OCM. 

ITV System 


In the ITV System a TOW mount is driven through rate command dynamics 
by the human using a handlebar controller. The system dynamics can be ap- 
proximated by the transfer function 


T(s) 


K 

s(Ts ^ 1) 


T = - sec, 
TT 


K= .1 


Since the handlebar is spring-loaded we assume that the motor-noise scales 
with commanded control input. The motor-noise parameters pertinent to ITV 
are 


V° “ [.OA, .011. p •« .03, p„. » .012 
u u EA 

In additiou it was found that a weighing Q. = .5 resulted in a best match 
between Model and data PSD ensembles. 

TABLE 3 - COMPARISON OF TEMPORAL TRACKING ERRORS, ITV SYSTEM 


Azimuth SD Elevation SD 

Data Modal Data Model 


1C 2E 0.112(0.035) 0.114(0.027) 0.092(0.023) 0.091(0.017) 
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There was data from only one ensemble for ITV, corresponding to a cross- 
ing target at 2Km with 0^ ~ 1. The target was movirig towards the gunner on 
a 40° angle. Table 3 gives the model-data comparisons for RMS tracking error 
in this stationary ensemble. A subjective comparison of a typical model-vs- 
data time history is given in Fig 10. The comparison of PSD ensemble sta- 
tistic of model and data is provided in Figs 11-12. Again, we find excellent 
agreement between OCM results and the field-test data. Note that this agree- 
ment is excellent not only for the PSD mean statistics, but also for the 2nd- 
order statistics that give an indication of the lun-to-run variability. 


CONCLUSIONS 

It has been demonstrated that the Optimal Control Model can be used to 
generate accurate facsimiles of target tracking error in various different 
anti-tank systems. While these results were not entirely unexpected, based 
on previous applications of the OCM, they are quite Interesting in that com- 
parisons have been made across several dimensions. By using the model to 
generate an ensemble, data and model ensembles can be studied, averaged and 
manipulated in similar manners, yielding slmllai results (at least to 2nd 
Older statistics) . 

The types of systems studied were quite varied, especially with regard 
to their manipulator characteristics. Thus, it was necessary to adjust the 
motor-noise paramccers in the OCM among systems. 

Further application of the OCM to anti-tank tracking systems is expected 
to refine these results, focus on dynamic inter-axis attentional allocation, 
and refine techniques for parameter value identification. 


ACKNOWLEDGEMENT 


The author would like to acknowledge the personnel at the Army Material 
Systems Analysis Agency (AMSAA), in particular, Mr. Patrick O'Neil and 
Mr. Patrick Corcoran, for their support and assistance in supplying what 
seemed like endless numbers of tracking runs for the various systems studied. 


REFERENCES 


1. Kleinman, D.L., "A Model for Generating Human Tracking Error in Anti- 
tank Systems", AI.PHATECH, Inc., Burlington, Mass . Interim Report on 
Contract DAAK11-80-C-0050, July 1981. 

2. Kleinman, D.L., "Monte-Carlo Simulation of Human Operator Response?', 
Technical Report TR77-1, Dept, of EECS, University of Connecticut, 1977. 


- 386 - 



! 



o o d 

(5ss/psj)’*n 





- 387 - 


il 













Ensenible 








Ensemble Mean 























